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Abstract 

We discuss the case when a fermion occupies an excited non-zero frequency level 
in the field of domain wall. We demonstrate that a solution exists for the coupling 
constant in the limited interval 1 < g < g ma x ~ 1-65. We show that indeed there are 
different branches of stable solution for g in this interval. The first one corresponds 
to a fermion located on the domain wall (1 < g < v / 2vr)- The second branch, 
which belongs to the interval \P1tt < g < gmax, describes a polarized fermion off the 
domain wall. The third branch with 1 < g < g ma x describes an excited antifermion 
in the field of the domain wall. 



1 Introduction 

In our previous paper [1] we studied the problem "domain wall + excited fermion". Initially 
the problem of the spectrum of a fermion coupled to the field of a static kink was discussed 
in Refs. J2JEHHE1E]- These papers were devoted mainly to a zero-frequency fermion bound 
by the domain wall. 
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Fermionic bound states in the field of external kink were studied in Ref. [7] for the 
case of the A0 4 -model and in Ref. [8] for the case of the sine-Gordon model. However, 
the authors of these publications have considered kink as given external field. As we shall 
demonstrate in this our work, the presence of the fermion changes drastically the kink 
profile. So the excitation spectrum for the problem "fermion coupled to kink" looks quite 
different from that calculated in the external field approximation. 

We studied the system of the interacting scalar (0) and fermion fields in two- 
dimensional space-time (1 + 1). In terms of dimensionless fields, coupling constant g and 
space-time variables (x,t), the Lagrangian density was taken in the form 



The equation of motion for scalar field 0(x, t) in the presence of a fermionic field \1/ reads: 



where ^ = (3 is the Pauli matrix, see (II. 5p . 

If the coupling of the scalar field to fermions is switched off, g = 0, the equation of 
motion ( II. 2p has a static solution called "kink", 



In three space dimensions this solution corresponds to a domain wall that separates two 
space regions with different vacua <f>± = ±1, see, e.g. (HI ED] for more details. 

Let us discuss the fermionic sector of the theory. After the substitution ^(x, t) = 
e~ t£t ip £ (x) the Dirac equation for the massless case reads: 




(1.1) 



d^4> - 20 + 20 3 = -gW, 



(1.2) 



4>k(x) = tanhx. 



(1.3) 




(1.4) 



where a x and (3 are the Pauli matrices, 




(1.5) 
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In equation (jl.4p . 



if> s (x) 




is the two-component spinor wave function. In terms of functions u £ (x) and v e (x) 
Eq. (jl.4p takes the form 

(1.6) 



^ + g<p(x)u e 
ax 



dv £ 

dx 



+ 



Substituting <p(x) = 4>k(x) = tanhx we finally get: 



d 2 



g{g 



dx 2 
d 2 v F 



dx 1 



cosh 2 x 6 



cosh x 



2 2\ 

9 K, 



(1.7) 



(^ 2 - *V 



This system describes the spectrum and eigenfunctions of the fermion in the external 
scalar field 4>k{x) = tanhx. Solutions of Eq. (jl.7p with e 2 > g 2 belong to the continuum 
and those with e 2 < g 2 to the bound states. The best known discrete mode is the so-called 
zero-mode solution (it is time-independent, e = 0): 



r(g+i/2) 



i 



cosh 9 x 



V 



o 



/ 



For zero-mode solution ( II. 8p the r.h.s. of Eq. (jl.2p ^>^> = 2uv = 0, so we conclude that 
the solution of the full problem in the form "0#(x) + zero-mode bound fermion ( II. 8p " is 
self-consistent. 

However, if the fermion occupies a level with e ^ 0, the r.h.s. of Eq. (II. 2p is different 
from zero. Hence the kink's profile has to be modified to fulfil Eq. (II. 2p . In our previous 
paper [1J we found one example of analytic solution for the excited fermion on a distorted 



domain wall, which indeed is self-consistent. 
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The plan of this our paper is the following. In Section 2 we study solutions for the first 
excited mode. We develop a simple variational procedure that allows to get a reasonable 
approximation to the solution for coupling constant g in the interval 1 < g < g m ax ~ 1-65. 
For gi = 2^2(2 — \^3) w 1.46 the approximation coincides with the exact solution, found 
earlier in [lj. In Section 3 we present detailed analysis for the obtained solutions. In 
particular, we demonstrate that, depending on the coupling constant, we get two solution 
branches for the excited fermion on domain wall. For the coupling constant in the interval 

1 < g < \/2k the excited fermion is practically localized on a distorted domain wall. In 
the interval of couplings \/2ii < g < g max the initial profile of the domain wall, Eq. (11.31) . 
is almost restored. We also found solution for an excited antifermion (i.e. solution with 
£ < 0) on the domain wall. This solution looks like a smooth function without cuts or 
jumps for any coupling constants g in the interval g G (1; g ma x]- 

At the same time for interval of the coupling \/2 < g < g ma x the two components of 
fermionic wave function behave very differently with respect to localization. Namely, the 
upper component u £ (x) of the spinor if) e (x) is located outside the domain wall, while the 
lower component v e (x) of the fermion wave function sits on the domain wall. So in this 
interval of g the domain wall separates in space the localizations of the upper and the 
lower components of the spinor. 

A general discussion and the summary of the results are presented in Section 4. 

2 Self-consistent solutions for fermion in the field of 
kink 

Let us look for a solution of the Dirac equation for the fermion in the field of a distorted 
kink (px — tanhax, where a is unknown real parameter to be determined from a self- 
consistency condition, as it will be discussed below. Introducing a new variable y = ax 
and a new parameter s = g/a, we get a system of equations for the fermionic wave 
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function, which formally coincides with that given by Eq. (II. 7ft : 



d 2 u E t s(s + 1) 
dy 2 cosh 2 y 

d 2 v £ i s(s — 1 



Uc-i = (e' 2 - s 2 W, 



(2.1) 



dy 2 cosh 2 y 6 ^ ^ 6 



Here e' — e/a. The wave function of the fermion for the first excited state in the field of 
the distorted kink (ft K reads (1 < s < +oo): 




aT(s - 1/2) 



/ , tanhax \ 

V2s - 1 — — ^ 

cosh ax 

1 



(2.2) 



cosh s_1 ax 

Substituting wave function f)2.2p and 4>k{%) — tanhax into equation (II. 2p . we get that 
for s = 2 and e' = \/3 the self-consistency equation (II. 2p is fulfilled if only the slope a 
satisfies the condition [T]: 



2a 2 -2 



-V3a 2 



a' 



2(2- V3). 



(2.3) 



So for the special case of s = 2 we obtain an analytic self-consistent solution for the 
problem of an excited fermion on the domain walj^. This solution corresponds to the 
value of coupling constant 



ff = 0i = 2V2(2->/3)«1.46. 



For arbitrary s G (1; +oo), s ^ 2 there is no analytic solution of equation (II. 2p in the 
form 4>(x) = tanhax. Substituting (12.21) into (11.21) we obtain the following constraint: 

2a 2 -2 = _^r (a - 1/2)^31 1 

V^Fr(s-l) cosh 2s - 4 ax 1 ; 

This constraint becomes an algebraic equation for a in the limit s — > 1, (and g = 1) with 



a = 1, and also for s = 2 (and g = gi) with a = a{2) = y 2(2 — v^)- The latter case 
corresponds to the exact solution discussed above, see Eq. (12. 3p . 



^^For s = 2 and e' = — v3 we get analytic solution for an antifermion on the domain wall. Its wave 
function is also given by Eq. 




Figure 1: slopes a(s) and /3(s). 

However, Eq. (12.41) may be used to determine the slope a(s) at x = for the fermionic 
wave function at arbitrary s G (1; +oo). The function a(s) is shown in Fig. 1. Note that 
for large s ^> 1 

ais) m 

s 

and lim a(s) = 0. 

s— s>+oo 

Inserting fermionic wave function (12. 2p with the so obtained a(s) into the r.h.s. of 
Eq. (II. 2p . we get the following equation for the scalar field <fi{x): 

d 2 <P , ^^_ a\s)sV2^lY{s-l/2) tanh(q( a ) x) 

dx 2+ Z( P W ~ 2v ^r( S -l) cosh 2s - 2 (a(s)xY [ ] 

This equation has no solutions for boundary conditions S (O) = 0, - = a(s), <fi — > 1 

with x — > +oo at arbitrary s G (1; +oo). The only exceptional cases are those with s = 1 
and s = 2, for which we found exact solutions discussed above. However, one may try to 
solve the equation of motion ( 12. 5 1) with modified boundary conditions for the field <fi s {x, 0): 

0.(0) = 0, **M = l3(s), Jim>(x) = l. (2.6) 



x— >+oo 
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Figure 2: profiles of scalar field 4> s (x) at different values of parameter s. Solid line - 
numerical solution of eq. (12. 5ft ; dashed line - approximation of the exact solution by 
4> s {%) — tai\h((3(s)x). Dotted line - the function 4>™(x) = tanh(a(s):r). 

We solved Eq. (12. 5p with boundary conditions (12. 6|) numerically, using the shooting 
method for the whole interval of parameter s G (l;+oo). The resulting function (3(s) 
is shown in Fig. 1. Note that (3(s) = a(s) at two points s — 1 and s — 2, for which our 
procedure reproduces exact solutions. In the region of small s G (1; 2) functions /3(s) and 
a(s) do not differ drastically. On the other hand, lim /3(s) = 1, what means that our 
solution for scalar field approaches undistorted kink in the limit of large s. 

Numerical solutions for the scalar field <p s (x) at different values of parameter s are 
given in Fig. 2, where we also show the profile of the function <p s {x) = tanh(/3(s)x). This 
figure illustrates that for all s G (1; +oo) the exact numerical solution is very close to 
tanh(/3(s)x). So one can say that with high accuracy the field configuration <p s {x) is a 
distorted kink with the slope /3(s) at the origin. 
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Figure 3: a) the energy of the fermion field e(s); b) total energy of the system "kink + 
excited bound fermion" E + (s) (solid line), the boundary line of the region "domain wall 
+ fermion in continuum" (dashed line). 

In terms of fields {(f) s (x); ^ s (x)} the procedure we used here is indeed a variational- 
type one. This approximate procedure reproduces exact solutions for s — 1, 2 and gives 
reasonable results in the whole range of parameter s. 

3 Stable branches of solution for a fermion coupled to 
domain wall 

The energy of a massless fermion in the field of the distorted kink <f> K (ct(s)x) is e(s) = 
\/2s — 1 a(s) and —e(s) for massless antifermion. The function e(s) is shown in Fig. 3a. 
Note that for large s ^> 1 e(s) goes to zero, e(s) ~ — 

The energy of the system "kink + bound fermion (antifermion)" may be well approx- 
imated by the following simple expression: 



The function E + (s) is shown in Fig. 3b. Notice that for each s G (1; +oo) the energy of 
the field configuration E + (s) is smaller than the energy of the configuration "undistorted 




(3.1) 
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Figure 4: coupling constant <? depending on parameter s. 

domain wall + fermion at rest in continuum". The energy of the latter configuration 
4 

Ef ree (s) = — + g(s) is also shown in Fig. 3b. Looking at this figure we conclude that the 
o 

field configuration "kink + excited fermion" is truly bound. 

Now let us consider the dependences of observables in terms of coupling constant 
g = g(s) = s a(s). The curve g(s) is shown in Fig. 4. As it is seen, the function g(s) 
reaches its maximum g rnax ~ 1.65 at s rn w 5.58 and falls down for large s ^> 1 reaching 
g as = \^2n 1.58 at s — >■ +oo. The asymptotic behaviour of g(s) is: 
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J( ,) = ^(i + - 

The crossing point g(si) = g as is also marked in Fig. 4 with si ~ 2.78. 

The energy of the system E±{g) as a function of coupling constant is shown in Fig. 5. 
Note that in the interval 1 < g < g as (1 < s < Si) the functions E + (g) and E^(g) 
are single-valued, whereas in the interval g as < g < g rnax {s\ < s < +oo) they are 
double- valued. So for each g G (g as ; g ma x] there are different branches of the solution for 
both E + (s) and E-(s). The first branch of function E + (s) corresponds to the interval 
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Figure 5: total energy E± of the system as a function of coupling constant g; stable 
branches of the solution are shown by solid lines; unstable branches are indicated by 
dotted lines. Horizontal dashed line demonstrates the energy of the udistorted kink. 

s £ [si; s m ]. The second branch corresponds to s G [s m ; +oo). As it is seen from Fig. 5, 
the first branch of solution E + (s) has larger energy E(g(s)) than the second one. That 
is why the first branch is unstable with respect to the decay into the second branch. The 
excess of energy could be emitted in the form of waves of scalar field (f). The unstable 
branch of the solution is drawn in Fig. 5 by dotted line. Analogous situation takes place 
for function E_(s). Stable branch of solution for function E_(s) is also shown in Fig. 5 
by solid line. 

The lower stable branch for function E + (s) corresponds to large s G [s m ; +oo). As it 
follows from Fig. 3a, in this region e(s) is small and e(s) — > at s — > +oo (or g — > g as ). 
That is why it is quite natural to call this branch the "quasi zero-mode" solution. 

So we got two distinct stable branches for field configurations: domain wall carrying 
an excited fermion. Let us discuss what these field configurations look like. 

The r.h.s. of Eq. (12.51) E(x, s) is proportional to the product u e '(x)v £ r(x). For different 
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F(x,s) 



0,40-i 




Figure 6: the right-hand side of eq. (12. 5p as a function of x for some values of s. 



values of parameter s the function F(x, s) is shown in Fig. 6. The maximum position for 

(s > 1), and 



each curve is at x T 



in 



Fornax: s) 



So in the limit of large s, i.e. for g — > g as = v2tt, the r.h.s. of Eq. (12. 5p becomes small 
and as a result the profile of the domain wall looks similar to the unperturbed kink 
4>k — tanhx. 

Let us clarify what the wave function of fermion looks like in the limit of large s ^> 1 . 
The position of the maximum of u 2 £ ,(x) is at x max = , , and 



u 2 £ i(x max ) — \ . 

V 7re 4 y/s 

Some profiles of u^,(x) are presented in Fig. 7a. We conclude that in the limit of large s 
the upper component of the fermionic wave function is located outside the domain wall. 
On the contrary, the maximum of the lower component v £ >(x) is at x = 0, i.e. on the 
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Figure 7: the functions u 2 E i{x) (figure a) and v E ,(x) (figure b) for some values of s. 
domain wall, for all values of s G (1; +00) and 




In the limit of large s, v E >(x) ~ exp(—g as \x\) (\x\ — > +00). The function v e /(x) behaves 
as a zero-mode solution in the limit of large s. Some profiles of v E ,(x) are presented in 
Fig. 7b. 

Thus, in the limit of large s the wave function of a fermion in the first excited state 
practically splits into two parts. The upper component of the wave function, u E i(x), which 
describes fermions with spin projection +1/2, is located outside the domain wall. The 
lower component of the wave function, v £ r(x), with spin projection —1/2, prefers smaller 
distances from domain wall. As a result, in the limit of large s the domain wall separates 
the maxima of the upper u e > (x) and the lower v £ > (x) components of the fermionic wave 
function in space. 

4 Conclusion 

We used a simple variational method and found a solution of the problem "domain wall 
+ excited fermion". We performed numerical simulations. For each coupling constant 
g G (l;<7max] the solution is a bound fermion with wave function (12. 2p in the field of 
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the distorted kink <pK = tanh(/3(s)x). The functions a(s) and (3(s) are shown in Fig. 1. 
Parameter s varies in the interval 1 < s < +00. This corresponds to a variation of 
coupling constant g in the limits 

max ~ 

1.65. 

No solution for a fermion in the first excited state on the domain wall was found for any 

9 9 max- 

For g G (l;g max ] we found two separate stable solution branches for function E + (s). 
For s G (1; Si), where Si = 2.78 and g(si) = g as = \^2tt 1.58, the fermion wave function 



is located on the distorted domain wall. For s = 2 (g = g\ = 2y2(2 — y/3)), belonging 
to this interval, our variational procedure reproduces the exact self-consistent solution of 
the problem, obtained in our previous publication [1]. 

For s G [si, s m ], where s m = 5.58, the solution is unstable. 

For the interval s G [s m ; +00), which corresponds to the variation of coupling constant 
g G [gmax'-, gas), the fermionic wave function consists of two different pieces approximately 

separated in space. The upper component u £ >(x) is located outside the domain wall and 

+00 

the lower one v e >(x) prefers smaller distances. Note that J \v £ i(x)\ 2 dx = — , which means 

—00 

that only one half of the fermionic wave function is located near the domain wall in this 
regime of large s. The profile of the domain wall for these large s, almost restores its 
original profile (fix = tanhx, as /3(s) — > 1 at s — > +00. 

Note that the energy of the fermion, e(s), is getting very small in the limit of large s 
(see also Fig. 3a). We conclude that in this limit, s > 1, the spectrum of the problem 
"domain wall + fermion" has two nearly degenerate states: the zero-mode ground state 

\/87T 

(e = 0) and the state with an excited fermion. The energy of the latter is s(s) ~ — j=r ^ 1, 

V s 

so this level may be classified as a quasi zero-mode solution. For g > g max only zero-mode 
solution survives. The picture is slightly different for an antifermion on the domain wall 
(function E-(s)). For this case the stable branch of solution corresponds to s G (1; s m ) 
(9 e (l;g max )). 
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Speaking of the practical applications of the results discussed above, we should first of 
all mention the problems of cosmology. In this connection we refer to the classical studies 
of Refs. PHH2UI3]. 

Note that a similar problem "kink + charged scalar field" at classical level was solved 
in Ref. |14] . A related problem of the interaction of a domain wall with a skyrmion was 
studied in Ref. [T5] . 

Some other properties of fermions in the fields of topological objects such as solitons 
and instantons one may also find in the monography by R. Rajaraman [16J. 
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